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IC.Sumnary/Notes LANDSAT (MSS) -imagery, aerial photography and the IMAGE- 
100 system were used in this study for demograpTrie estimations and for 
monitoring of urban growth. Two sets of urban test sites, one with 35 
oities and one with 70 cities Were selected in the State of Sao Paulo. A\ 
high degree of aolinearity (0.96) was fovmd between urban areal 
measurements taken from aerial photographs and LANDSAT (MSS) imagery. 
High coefficients were observed when census data was regressed agcdnst 
aerial information (0.95) and LANDSAT data (0.92). The validity of 
population estimations was tested by regressing three urban varieties, 
against three classes of cdties. The results supported the effectiveness 
of LANDSAT to estimate large city papulations with diminishing 
effectiveness as urban areas decrease in size. The IMAGE-100 system 
computed urban growth between October, 1972 and June, 1976. Residential 
and industrial sectors in Sdo Jose dos Campos showed growth along the 
Sdo Paulo to Rio de Janeiro urban "corridor" . Further research will 
refine the LANDSAT system's precision to analyse and predict 
urban growth in the State of Sdo Paulo, although the methodology should 
be of use to any other region, with adequate adaptation. 

- ■ . ■■ ' - 

OpmarKs ^bis paper will be presented at the "Simposio Intemaeional ^ 
varcevcion Bemota aplieada a Demografia y Uso actual de la Tierra", in 
topaz, Bolivia, (28-30 November). 



1 - INTRODUCTION 


The utilization of LANDSAT data for demographic research 
and for the monitoring of urban growth is one of the most challenging 
aspects of applications in Remote Sensing. 

The Brazilian government conducts a demographic census 
every 10 years to monitor population growth. However, inter-census 
estimations are considered to be inaccurate owing to the rapid growth of 
urban areas. 


The purpose of this study is to develop a viable 
methodology using LANDSAT data as a primary source for population 
estimation and assessment of urban growth in d:he intercensal period. 

This research in the use of low resolution imagery is an 
attempt to verify existing relationships between spectral reflectance 
characteristics and variations in urban population. 

The specific objectives of this study are to: 1) identify 
and quantify the size of individual urban areas using information 
obtained from LANDSAT and conventional imagery, and the IMAGE-100 system; 
2) make comparisons between the results obtained from LANDSAT data 
(scale = 1:500,000} and those obtained from conventional aerial 
photographs (scale = 1:25,000); 3) evaluate LANDSAT data as a monitor 
of the expansion of urban growth, and 4) analyze the relationship 
between population data, based on the "Demographic Census" of 1970, and 
urban area measurements obtained from LANDSAT positive imagery in 
spectral channel 5. 

The potential of LANDSAT data for population estimates 
has been examined in previous papers of Forestland Mendonga (1974). 

Other working papers on the subject include: Ellefsen et al (1974), 

Hsu (1973), Durland (1974) and Durland (1975), Reining (1974), 

Reining and Egbert (1975), Christenson et al (1975) and others. 
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1.1 - DEMOGRAPHIC FEATURES OF THE STUDY AREA 

The State of Sao Paulo was selected for this research 
for reasons of high urban growth and its geographic situation as a vital 
region of development and organization in Brazil (Figure 1). 



Fig. 1 - State of Sao Paulo: Study Area 
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The population of the State of Sao Paulo encompasses 
about 19..0/i of the total Brazilian population according to the national 
"Demographic Census” of 1970. This figure is impressive considering 
that the state's area (247,898 km^) represents only 2.95^ of the 
country's total area (8,511,955 km^). The demographic growth of the 
State of Sao Paulo is more a result of incoming migrants than of natural 
population growth. Moreover, the industrial development of the Sao Paulo 
metropolitan area and of other regional urban centers are important 
factors in the urban development of the whole State, 

Porcaro and Oliveira (1976) report that urbanization is 
a function of population concentration in urban centers. Table 1 shorn 
urban and total population growth from 1940 - 1950 and from 1960 - 1970 
for the State of Sao Paulo and for Brazil. 


TABLE 1 

URBAN AND TOTAL POPULATION GROWTH 



URBAN POPULATION (t) 

TOTAL POPULATION («) 


1940-1950 

1950-1960 

1960-1970 

1940-1950 

1950-1960 

1960-1970 

state of Sao Paulo 

51.6 

69 5 

80 3 

27 2 

42 0 

37.0 

Brazil 

45 8 

70.3 

62,9 

25.0 

36 7 

31.3 


Porcaro and Oliveira (1976,) concluded that the ratio between urban 
population and total population for a region is a primary indicator of 
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the degree of urbanization. Table '2 shows urbanization rates for the 
State of Sao Paulo and for Brazil from 1940 to 1970, based on 
“Demographic Census" data. 


TABLE 2 

URBANIZATION RATES 



1940 

1950 

1960 

1970 

State of Sao Paulo 

44. n 

52.8^ 

62.8% 

80.3% 

Brazil 

31 .2% 

36.1% 

45.8% 

55.9% 


2 - MATERIALS AND METHODS 
2.1 - MATERIALS 


This study was carried out using LANDSAT system products. 
The description and specification of these products can be found in the 
NASA Data User Handbook (1972). 

A visual interpretation of urban areas was conducted 
using black and white LANDSAT positives at scales of 1:1,000,000 and 
1:500,000. 


The tonal response of channel 5 (0.6 vm - 0.7 um) was 
found to be the most suitable for visual analysis of urban areas. The 
tonal contrast between urban and surrounding areas, observed in this 
channel, provided an excellent base for the delineation of urban 
boundaries. CCT were employed for automatic interpretation using 
the IMAGE-100 multispectral analyzer (GE, 1975). Black and white, 
conventional aerial photographs (scale = 1:25,000) were used to obtain 
areal information for comparison with LANDSAT data. Topographic charts 
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(1:1,000,000), transportation maps (1:1,000,000) and a packaged computer 
program were used as supplementary materials. 

2.2 - METHODS 


Two sets of urban study areas were identified on LANDSAT 
imagery in accordance with Carter's and Stone's methodology (1974). 
Anderson's et al (1.972) criteria were applied to delineate urban areas 
on aerial photographs and LANDSAT imagery. Included within the urban 
limits are industrial, commercial, institutional, and recreational land- 
~use types, service areas and transportation networks. Surface transport 
links among cities were collected from the imagery and existing maps, 
and, tested according to Ogrosky's methodology (1975). 

Urban population data, for the years 1973, 1975 and 1976 
were based upon “Demographic Census" for 1970 using the mathematical 
formula of geometrical growth: 

P(t)= P (0) (l + i)"' 


where: 


P (0) = lasted census data existing for the number of inhabitants; 
P (t) == population at year t (t = year of the projection); 
i = annual growth rate; 
n = number of years between P (t) and P (0). 

Simple and multiple Linear Regressions were applied to 
the data for: 1) data obtained from 70 cities, and, 2) data separated 
into three classes based upon population size. Class limits were based 
upon Sigueira’s methodology (1976) and adapted for the State of Sao 
Paulo, These are: a small city class composed of 35 cities with 

populations less than 30,000 inhabitants; a medium city class with 25 
cities and less than 100,000 inhabitants, and, a large city class with 
10 cities and populations greater than 100,000 inhabitants. 
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LAMDSAT data processed from the IMAGE- 100, by: 1) spectral 
signature study of the urban areas, 2) single cell classification, and 
3) monitoring of urban growth in two different periods. 

3 - RESULTS AND DISCUSSION 


Table 3 shows urban population data and urban area 
measurements obtained from aerial photographs and LANDSAT imagery for 
35 cities. Areal data from LANDSAT are usually overestimated when 
compared with data obtained from conventional photographs. The 
correlation coefficient of (0.96) obtained between aerial and LANDSAT 
imagery shows the compatibility of these data sources. The correlation 
coefficient between census data and aereal data from conventional 
photographs is (0.95). This value, which is slightly greater than the 
coefficient of (0.92) obtained for the correlation between census data 
and urban LANDSAT imagery, supports the suitability of satellite data 
as a source for estimating populations. 

Population estimates were based on an enlarged sample of 
70 cities provided by recent LANDSAT imagery (1975/1976). Figures 2, 3, 
and 4 show urban areas identified on three LANDSAT Images located in the 
central east and southeast region of the State of Sao Paulo. Only well 
defined urban areas could be used as test sites. Figure 5 shows the 
location of the test sites. The numeration for each city is included in 
Appendix A which presents urban population (P), urban area (A) and 
surface transport (L) data. Table 4 shows the results for the simple 
linear regression and Figures 6, 7, 8, and 9 illustrate dispersion 
diagrams for the four regression runs presented in Table 4. The results 
show that the best fit was obtained for the "total sample" with a r^ 
coefficient of 0.9418. The division of the sample into three classes did 
not improve the accuracy of the population estimation model. A test of 
the correlation coefficients was applied to the values among the three 
classes of cities. The results show no significant differences among the 
coefficient values at a p = 0,01 level. A multiple linear regression was 
run, where the surface transportation linkage index was introduced 
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T ABLE 3 

POPULATION AND URBAN AREA DATA (1973) FOR 35 TEST SITES 


KO 

CITIES 

POPULAFION 

URBAN AREA 
AERIAL PHOTOGRAPH 
(kra^) 

URBAN AREA 
LANDSAT 
(km^) 

1 

itirapiMa 

4,556 

0.76 


2 

AVANHANOAVA 

3,997 

0.81 


3 

VINHEDO 

7,944 

1.32 

2.16 

4 

BUPITAMA 

6,550 

1.53 

1.60 

5 

ITATIBA 

22,474 

1.90 

4.16 

6 

SUHARE 

13,751 

2.30 

3 24 

7 

miranoopolis 

11 ,677 

2.45 

2.52 

8 

SANTO ANASTRCIO 

14,144 

2 70 

3 40 

9 

SRO ROQUE 

17,178 

2.85 

2.90 

10 

SALTO 

20,264 

2.90 

3.00 

11 

ITAPIRA 

26,959 

3.15 

3.48 

12 

PEREIRA BARRETO 

20,364 

3.37 

•5 00 

13 

MOCOCA 

22,347 

3 64 

4 00 

14 

Sta BARBARA D'OESTE 

24,127 

3.92 

3.44 

15 

VALINIIOS 

22,083 

3.92 

4.60 

16 

MOCI GUASSU 

34.P20 

4.04 

6.44 

17 

LEME 

26,132 

5. 11 

5.92 

18 

BRAGANgA PAULI STA 

40,472 

5.26 

4.72 

19 

ARARAS 

44,478 

5.29 

7.96 

20 

BIRI6UI 

28,698 

5 37 

7.64 

21 

PENRPOLIS 

26,661 

5.65 

7.00 

22 

PIRASSUNUNGA 

27,046 

5.67 

5.92 

23 

JAU 

42,678 

5.98 

4.56 

24 

ITU 

48,959 

6.40 

12 60 

25 

AMORADINA 

48,675 

7.06 

9.92 

26 

BOTUCATU 

55,835 

9.50 

12.60 

27 

RIO CLARO 

72,819 

13.81 

14.80 

28 

PRESIDCNTC PRUDCNTE 

98,696 

14.80 

14.84 

29 

PIRACJCABA 

134,217 

15.09 

17.40 

30 

MOGI OAS CRU2LS 

95,183 

15.94 

17 03 

31 

SRO CARLOS 

79,565 

17.21 

17.80 

32 

ARAQATUBA 

91,890 

18.56 

18.76 

33 

ARARAQUARA 

87,045 

19 45 

19.88 

34 

SOROCABA 

176,742 

22.40 

28.70 

35 

RIUEIRRO PRCTO 

205,597 

26.45 

21.24 











priolfbiij 


Bototait 


MINAS 


' Virodourv 


Jordkwpolis 


,XPo<lt«l 


Bebcdoi 


■Ri^eff 

Pr«»o 


BorrinlM 


J doboticflbfl^ 

jC&ndMo^ddrisot* 

■■■' A_*nr-f^'iW 

^ TnrkMirat^naf^ md 


MRibtirJk' 

Benito 


Concho! 


Sonto <i«rtrw4«m 


PROCeSSAEO 24N0V75 


IS 


10MAI75 C S21-32^W047-49 N S2l-32^W047-44 MSS 5 SOL EL33 A2049 190 14234 P L2 3 01 IHPE LAHDSAT- 175130- 121903-5 


LEGEND 

CITIES 

RESERVOIRS 

INTERSTATE BOUNDARIES 

RIVERS 

HIGHWAYS 


07 POOR 





- 10 - 



tt 

Fui'U 


MINAd 


Cochoelra 


LOfftlO 


fkptttdo 


Pin;|o0'onh4h9abo ^ 


Tauh^ii 


^^Quori 





25JUN76 ORB 150 PT 28 Cl S23-04/U045-23 5SS 5 

Nl S23-02/U045-14 


SOLI EL22 B2049 B 190 LC 19979 INPE^LftNDSPT 176177-114204-5 

R02 N -BHRSIL- 28M0R77 CEN« 017 


LEGEND 


CITIES 

RESERVOIRS 

INTERSTATE BOUNDARIES 

RIVERS 

HIGHWAYS 


ORIGINAL PAGE IB 
OF POOR QUALITY 


Fio. 4 - Delimitation of Urban Areas of the Sao Jose dos famnos Reoion 






- 12 - 


TABLE 4 

SIMPLE LINEAR REGRESSION 



rt 

r* 

y = a + bx 

Sy.x (no inhabitants) 

Small Cities 

(less than 30,000 inhabitants) 

35 

0.7289 

Pi = 4,061 + 4,507 (Al) 

3,844 

Hediun Clues 

(30,000-100,000 inhabitants) 

25 

0Jb147 

Pi = 14,045 + 4,444 (Ai) 

8,676 

Urge Cities 

(more than 100,000 inhabitants) 

10 

0.8529 

Pi = 16,613 a 6,337 (Ai) 

46,630 

TOTAL SAMPLE 

70 

0.941S 

Pi = 1,186 +'6,672 (Al) 

13,355 


into the model. Holtz' et al (1969), Ogrosky (1975) and Foresti (1976) 
found this variable to be a useful estimator of population after the 
“urban area" variable. Table 5 shows the results of the regression run 
where the addition of the transportation index increased the coefficients 
of determination (R^): for sma.ll cities (0.0020), medium cities (0.0135), 
and, large cities (0.0026). Ogrosky (1975) reported an increment of 
0,0068 in the coefficient when the surface transportation index was 
included in the multiple algorithm, The standard error (S ) showed an 

* y t A 

increment for all classes of population size except for the medium class. 

Appendix B and C provide estimates of population and 
error estimates found for each of the 70 urban areas based upon simple 
Linear Regression and Multiple Linear Regression, respectively. 

It must be pointed out that a small error factor of about 
8 . 0 % is inherent in our projections of population growth for 1976 based 
on the Demographic Census of 1960/1970, 




POPULATION tlN THOUSANDS) 



POPULATION (IN THOUSANDS) 
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Fig. 7 - Relationship Between Population and Urban Area for 
Medium Cities 



POPULATION ( IN TROUSANDS) 
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fig 8 - Relationship Getv/een Population and Urban Area for 
Large Cities 


POPULATION (IN THOUSANDS} 
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TABLE 5 

MULTIPLE LINEAR REGRESSION 



K 

■i 


Sy.x 

(no. inhabitants) 

Small Cities 

(less than 3,000 inhabitants) 

35 

0 7309 

Pi = 3.224 + 4,462(Ai) + 63(Li) 

3,889 

Medium Cities 

(30,000 - 100.000 inhabitants) 

25 

0 8282 

Pi =25.481 <4,306(Al} - 625(Li) 

8,542 

Large Cities 

(more than 100,000 inhabitants) 

10 

0 85SS 

Pi = 10,495 4 6,952 (Ai) +1,217(Li 

4.907 

TOWL SAMPLE 

70 

0 9413 

Pi =-2.231 1 6.638(Ai) + 66(Li) 

18/88 


Urban area data collected from LANDSAT imagery also 
presented problems. Among the difficulties encountered were such factors 
as shape (linear or star shaped areas). Also difficult to analyze are 
areas located in regions of rough topography where shadow effects are 
created. Cities surrounded by crop belts tend to obscure tonal features 
unless coverage is provided for the wet season, Murray (1975). In 
addition, cities located on or near the coast often exhibit overlapping 
reflectance between city and beach. 

The decrease of the error associated with the urban 
population can be lowered using imagery taken at the same period of the 
demographic census. 

Greater precision, for collecting urban area data, can 
be obtained using higher resolution images and by employing correction 
factors through magnification or enhancement techniques. 
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An urban expansion study of Sao Jose dos Campos, using 
the IMAGE-100 system was conducted by the following procedure: 

1} Examination of LANDSAT (MSS) black and white positives (scale = 

- 1:500,000); 

2) Visual monitoring of black and white CRT images in each of four 
channels (4, 5, 6, 7) at a scale of 1:60,000; 

3) Selection of "training sites" and development of an urban land 
use classification system; 

4) Examination of spectral signatures for class I land-use types 
(Anderson et al, 1972); 

5) Analysis of alphanumeric computer printouts and redefinitions of 
"training sites"; 

6) Utilization of an auxiliary program algorithm to compare urban 
growth between the periods of October 1972 and June 1976. 

The study of urban structure is among the most complex in 
Remote Sensing and requires the combination of several basically different 
ground cover classes into a single land use classification system 
(Odenyo and Pettry, 1977). Because city componentes are not uniformly 
distributed within any one area, they are spectrally multi-modal. Such is 
the case for Sao Jose dos Campos where the difference in spectral 
signatures, for each CRT image, are probably caused by seasonal 
differences in data acquisition. For example, June 1976 (dry season) data, 
presented higher reflectance in all channels compared to the signatures 
of other images examined. Residential and industrial areas presented mixed 
spectral responses which can be confused with the response from crop areas, 
characterized by bare soil and vegetation. This fact makes spectral 
discrimination between urban and rural areas difficult when automatic 
analysis is performed. Figure 10 shows the output from the interactive 
digital system in alphanumeric printout characters. The printout 
represents a data element or pixel which corresponds to the cell 
resolution of the satellite. Figure 10 represents the Sao Jose dos Campos 
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Fig. 10 - IMAGE-100 Generated Alphanumeric Printout for Sao Jose dos Campos 
Urban Area in 1972 and for Urban Growth up to 1976 


The results of this study provide for a viable approach 
to the problem of population estimation and urban expansion and has 
future implications in the area of computer mapping and land use 
assessment. 
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BASIC DATA FOR SELECTED TEST SITE 


LINKS 











- A. 2 - 


NQ 

TEST SITE 

POPULATION 

AREA 

(km^) 

LINKS 

29 

JABOTICABAL 

35,619 

5.04 

25 

30 

JACAREI 

69,665 

8.28 

19 

31 

JAU 

49,076 

5.80 

22 

32 

JUNDIAI 

215,568 

27,00 

30 

33 

LEME 

34,424 

6.20 

17 

34 

LIMEIRA 

108,550 

23.00 

32 

35 

LORENA 

51 ,750 

9.08 

20 

36 

MOCOCA 

26,111 

5,20 

17 

37 

MOGI DAS CRUZES 

112,856 

2Q.64 

24 

38 

MOGI GUAQU 

51 ,490 

9.12 

20 

39 

MOGI MIRIM 

37,773 

4.28 

20 

40 

MONTE ALTO 

18,585 

3.76 

14 

41 

PIEDADE 

8,859 

1.20 

12 

42 

PINDAMONHANGABA 

36,400 

8.36 

29 

43 

PIQUETE 

14,250 

3.36 

15 

44 

PIRASSUNUNGA 

31,976 

5.96 

24 

45 

PONTAL 

8,132 

1.16 

13 

46 

PORTO FELIZ 

15,409 

1.76 

13 

47 

PORTO FERREIRA 

20,576 

4.60 

26 

48 

RIBEIRSO PRETO 

253,122 

23.04 

! 

35 

49 

RIO CLARO 

83,657 

15.08 

27 

50 

RIO DAS PEDRAS 

13,220 

1.68 

20 

51 

SALTO 

24,906 

3.68 

18 

52 

SANTA BARBARA 

31 ,096 

6.52 

23 

53 

SANTA CRUZ OAS PAL 
MEIRAS 

11,322 

1.32 

17 - 

54 

SANTA GERTRUDES 

5,548 

1.32 

10 

55 

SAO CARLOS 

92,784 

19.60 

29 

56 

SAO JOSE DOS CAMPOS 

225,463 

30.40 

32 

57 

SAO ROQUE 

19,686 

2.96 

16 

58 

SAO SIMAO 

9,037 

1.24 

24 








- A. 3 - 


NO 

TEST SITE 

POPULATION 

AREA 

(km^) 

LINKS 

59 

SSO VICENTE 

156,384 

18,60 

20 

60 

SERRANA 

10,779 

1.52 

10 

61 

SERRA NEGRA 

7,573 

1.72 

10 

62 

SERTAOZINHO 

28,611 

5.16 

15 

63 

SOCORRO 

9,506 

1.36 

16 

64 

SOROCABA 

216,676 

32.20 

27 

65 

SUMARE 

23,059 

4.32 

14 


66 

TAQUARITINGA 

22,401 

3,36 

67 

TAUBATE 

129,163 

17.80 

68 

UBATUBA 

15,290 

3.32 

69 

VALINHOS 

31 ,387 

4.92 

70 

VOTORANTIN 

43,270 

3.48 









APPENDIX B 


RELATIONSHIP BETWEEN ACTUAL AND ESTIMATED VALUES 


USING SIMPLE LINEAR REGRESSION 


ACTUAL POPULATION 


ESTIMATED POPULATION 


DIFFERENCE 


96,271 

23,703 

30,153 

100,377 

57,878 

25.068 
9,841 

25,132 

37.069 
50,427 
32,184 
17,623 

485,409 

13,785 

19,777 

14,314 

18,422 

13,410 

6,068 

57,130 

28,954 

69,811 

13,211 

31 .070 
33,629 

4,188 

29,228 


104,767 
17,762 
20,164 
149,871 
60,998 
32,441 
8,955 
29,772 
36,978 
41 ,248 
34,576 
12,691 
496,557 
19,363 
25,502 
13,492 
16,961 
15,894 
9,222 
48,187 
15,360 
75,676 
8,955 
34,576 
33,508 
3,617 
33,508 








- B.2 - 


N9 

ACTUAL POPULATION 

ESTIMATED POPULATION 

DIFFERENCE 

28 

47,950 

56,994 

- 9,044 

29 

35,619 

32,441 

3,177 

30 

69,665 

54,059 

15,605 

31 

49,076 

35,512' 

11,563 

32 

215,568 

178,962 

36,605 

33 

34,424 

40,181 

" 5,757 

34 

108,550 

152,273 

- 43,723 

35 

51 ,570 

59,396 

- 7,646 

36 

26,111 

33,508 

- 7,397 

37 

112,856 

136,527 

- 23,671 

38 

51 ,490 

63,667 

- 12,177 

39 

37,773 

27,370 

10,402 

40. 

18,585 

23,900 

- 5,315 

41 

8,859 

6,820 

2,038 

42 

36,400 

54,592 

- 18,192 

43 

14,250 

21 ,232 

- 6,982 

44 

31 ,976 

38,579 

~ 6,603 

45 

8,132 

6,553 

1,578 

46- 

15,409 

10,556 

4,852 

47 

20,576 

29 ,505. 

- 8,929 

48 

253,122 

152,540' 

100,581 

49 

83,657 

99,429 

- 15,772 

50 

13,220 

10,022 

3,197 

51 

24,906 

23,367 

1,538 

52 

31 ,096 

42,316 

- 11,220 

53 

11,322 

7,620 

3,701 

54 

5,548 

7,620 

- 2,072 

55 

92,784 

129,588 

- 36,804 

56 

225,463 

201 ,647 

23,815 

57 

19,686 

18,563 

1 ,122 

58 

9,037 

7,087 

1,949 








- B.3 - 


NQ 

ACTUAL POPULATION 

ESTIMATED POPULATION 

DIFFERENCE 

59 

156,384 

122,915 

33,468 

60 

10,779 

8,955 

1,823 

61 

7,573 

10,289 

- 2,716 

62 

28,611- 

33,242 

- 4,631 

63 

9,506 

7,887 

1,618 

64 

216,676 

213,657 

3,018 

65 

23,059 

27,637 

- 4,578 

66 

22,401 

21 ,232 

1,168 

67 

129,163 

117,578 

11,584 

68 

15,290 

20,965 

- 5,675 

69 

31,387 

31 ,640 

253 

' 70 

43,270 

22,032 

21,237 







APPENDIX C 


RELATIONSHIP BETWEEN ACTUAL AND ESTIMATED VALUES 


USING MULTIPLE REGRESSION 






ACTUAL POPULATION 


96,271 

23,703 

30,153 

100,377 

57,878 

25.068 
9,841 

25,132 

37.069 
50,427 
32,184 
17,623 

485,409 
13,785 
19,777 
' T4',314 
18,422 
13,410 
■6,068 
57,130 
28,954 
69,811 
13,211 
31 ,070 
33,629 
4,188 
29,228 


ESTIMATED POPULATION 

DIFFERENCE 

105,162 

- 8,891 

17,677 

6,025 

20,199 

9,953 

150,102 

- 49,725 

61 ,022 

- 3,144 

31 ,885 

- 6,817 

8,783 

1,057 

29,494 

- 4,365 

37,785 

716 

41 ,505 

8,921 

34,735 

- 2,551 

12,500 

5,122 

496,009 

- 10,600 

19,138 

- 5,353 

25,179 

- 5,402 

13,825 

488 

16,682 

1,739 

16,809 

~ 3,399 

8,586 

- 2,518 

48,541 

8,588 

15,287 

13,666 

75,824 

- 6,013 

8,387 

4,823 

35,131 

- 4,061 

33,673 

44 

2,944 

1,243 

33,211 

- 3,983 













- C.2 - 


NO 

ACTUAL POPULATION 

ESTIMATED POPULATION 

28 

47,950 

57,105 

29 

35,619 

32 ,875 

30 

69,665 

53,986 

31 

49,076 

37,722 

32 

215,568 

178,978 

33 

34,424 

40,047 

34 

108,550 

40,047 

35 

51 ,750 

59,363 

36 

26,111 

33,409 

37 

112,856 

136,364 

38 

51,490 

63,611 

39 

37,773 

27,500 

40 

18,585 

23,652 

41 

8,859 

6,529 

42 

36,400 

55,178 

43 

14,250 

21 ,063 

44 

31 .976 

38,916 

45 

8,132 

6,327 

46 

15,409 

10,310 

47 

20,576 

30,020 

48 

253,122 

153,021 

49 

83,657 

99,654 

50 

13,220 

10,241 

51 

24,906 

23,385 

52 

31 ,096 

42,567 

53 

11,322 

7,653 

54 

5,548 

7,191 

55 

92,784 

129,790 

56 

225,463 

201 ,680 

57 

19,686 

18,474 

58 

9,037 

7,584 


DIFFERENCE 


- 9,155 
2,743 

15,678 
11.353 
36 ,589 

- 5,623 

- 5,623 

- 7,613 

- 7,298 

- 23,508 

- 12,121 

10,272 

- 5,067 
2,232 

- 18,778 

- 6,813 

- 6,940 
1,804 
5,098 

- 9,444 

100,100 

- 15.997 

2,978 
1 ,520 

- 11,471 

3,668 

- 1 ,643 

- 37,006 
23,782 

1,211 

1,452 








- C.3 - 


NO 

ACTUAL POPULATION 

ESTIMATED POPULATION 

DIFFERENCE 

59 

156,384 

122,558 

33,825 

60 

10,779 

8,519 

2,259 

61 

7,573 

9,846 

- 2,273 

62 

28,611 

33,011 

- 4,400 

63 

9,506 

7,853 

1 ,652 

64 

216,676 

213,299 

3,376 

65 

23,059 

27,369 

- 4,310 

66 

22,401 

21 ,129 

1,271 

67 

129,163 

118,238 

10,924 

68 

15,290 

20,995 

- 5,705 

69 

31 ,387 

31 ,665 

229 

70 

43,270 

21,793 

21 ,m 








